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Abstract. We derive the high frequency limit of the Helmholtz equation with source 
term when the source is the sum of two point sources. We study it in terms of Wigner 
measures (quadratic observables) . We prove that the Wigner measure associated with the 
solution satisfies a Liouville equation with, as source term, the sum of the source terms 
that would be created by each of the two point sources taken separately. The first step, 
and main difficulty, in our study is the obtention of uniform estimates on the solution. 
Then, from these bounds, we derive the source term in the Liouville equation together 
with the radiation condition at infinity satisfied by the Wigner measure. 

AMS subject classifications. 35Q60, 35J05, 81S30. 

1. Introduction 

In this article, we are interested in the analysis of the high frequency limit of the 
following Helmholtz equation 

Cv Til 

(1.1) -i—u £ + Au £ + ^-u £ = S £ (x), x€R 3 

with 

S £ (x) = S £ (x) + Sl<x) = iSo(J) + ^i(^) 
where qi is a point in R 3 different from the origin. 

In the sequel, we assume that the refraction index n is constant, n(x) = 1. 

The equation l|l.l|l modelizes the propagation of a source wave in a medium with 
scaled refraction index n(x) 2 /e 2 . There, the small positive parameter e is related to 
the frequency u> = 557 of u £ . In this paper, we study the high frequency limit, i.e. 
the asymptotics e — ► 0. We assume that the regularizing parameter a £ is positive, 
with a £ — ► as e — * 0. The positivity of a e ensures the existence and uniqueness 
of a solution u £ to the Helmholtz equation Ijl.lfl in L 2 (R 3 ) for any e > 0. 

The source term S £ models a source signal that is the sum of two source signals 
concentrating respectively close to the origin and close to the point q± at the scale 
s. The concentration profiles So and Si are given functions. Since e is also the 
scale of the oscillations dictated by the Helmholtz operator A + , resonant inter- 
actions can occur between these oscillations and the oscillations due to the sources 
Sq and Sf. On the other hand, since the two sources are concentrating close to 
two different points in R 3 , one can guess that they do not interact when e — > 0. 
These are the phenomena that the present paper aims at studying quantitatively. 
We refer to Section 3 for the precise assumptions we need on the sources. 

In some sense, the sign of the term —ia £ eu £ prescribes a radiation condition at 
infinity for u £ . One of the key difficulty in our problem is to follow this condition 
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in the limiting process s — > 0. 

We study the high frequency limit in terms of Wigner measures (or semi-classical 
measures). This is a mean to describe the propagation of quadratic quantities, like 
the local energy density |w e (x) | 2 , as e — > 0. The Wigner measure [i(x, £) is the 
energy carried by rays at the point x with frequency £. These measures were intro- 
duced by Wigner and then developed by P. Gerard [H| and P.-L. Lions and T. 
Paul H] (see also the surveys [3] and [7|). They are relevant when a typical length 
e is prescribed. They have already proven to be an efficient tool in the study of 
high frequencies, see for instance [5], 0] for Hclmholtz equations, P. Gerard, P. A. 
Markowich, N.J. Mauser, F. Poupaud 7 for periodic media, G. Papanicolaou, L. 
Ryzhik for a formal analysis of general wave equations, L. Erdos, H.T. Yau [H] 
for an approach linked to statistical physics, and L. Miller 1Q. for a study in the 
case with sharp interface. 

Such problems of high frequency limit of Helmholtz equations have been studied 
in Benamou, Castella, Katsaounis, Perthame 2 and Castella, Perthame, Run- 
borg [3] . In [2] , the authors considered the case of one point source and a general 
index of refraction whereas in [3], they treated the case of a source concentrating 
close to a general manifold with a constant refraction index. In the present paper, 
we borrow the methods used in both articles. 

In the case of one point source, for instance Sq only, with a constant index of re- 
fraction, it is proved in that the corresponding Wigner measure /io is the solution 
to the Liouvillc equation 

OVo(x, + £ • Vx/zofo = Qo(x, = -^5(x)5(\e ~ 1)|^(£)| 2 , 

the term + meaning that [i is the outgoing solution given by 

Vo(x,Q = Qo(x + t€,g)dt. 

J —oo 

In particular, the energy source created by Sq is supported at x = 0. Similarly, 
the energy source created by the source Sf is supported at x = q\. Thinking of 
the orthogonality property on Wigner measures, one can guess that the energy 
source generated by the sum Sq + Sf is the sum of the two energy sources created 
asymptotically by Sq and Sf. 

Indeed, we prove in this paper that the Wigner measure /i associated with the 
sequence (u £ ) satisfies 

(1.2) 0+^,0 +£-V*/i(z,0 = Qo + Qi, 

where Qo and Qi are the source terms obtained in J5] in the case of one point 
source. However, our proof does not rest on the mere orthogonality property. 

Let us now give some details about our proof. Our strategy is borrowed from [5] . 
First, we prove uniform estimates on the sequence of solutions (u e ). We also study 
the limiting behaviour fo the rescaled solutions e^~u e {ex) and e^~u £ (qi + ex). 
The obtention of these first two results is the key difficulty in our paper. It relies 
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on the study of the sequence (a £ ) such that 

— ia £ ea e + Aa £ + a £ = Si — — ^ . 

Using the explicit formula for the Fourier tranform of a s , we prove that a e is uni- 
formly bounded in a suitable space and that a £ — > as e — > weakly. We would 
like to point out that our analysis, based on a study in Fourier space, strongly rests 
on the assumption of a constant index of refraction. 

Second, our results on the Wigner measure then follow from the properties proved 
in [2]. They are essentially consequences of the uniform bounds on (u £ ): we write 
the equation satisfied by the Wigner transform associated with (u £ ), and pass to 
the limit e — > in the various terms that appear in this equation. The only difficult 
(and new) term to handle is the source term. 

Third, we prove an improved version of the radiation condition of 2 . Our argu- 
ment relies on the observation that fi is localized on the energy set {|£| 2 = 1}, a 
property that was not exploited in [2]. 

The paper is organized as follows. In Section 2, we recall some definitions and 
state our assumptions. Section 3 is devoted to the proof of uniform bounds on the 
sequence of solutions (u £ ) and of the convergence of the rescaled solutions. Then, 
in Section 4, we establish the transport equation satisfied by the Wigner measure 
fj, together with the radiation condition at infinity. In the appendix, we recall the 
proof of some results established in [2] that we use in our paper. 

2. NOTATIONS AND ASSUMPTIONS 

In this section, we recall the definitions of Wigner transforms and of the B, 
B* norms introduced by Agmon and Hormander £Q for the study of Helmholtz 
equations. Then, we give our assumptions. 

2.1. Wigner transform and Wigner measures. We use the following definition 
for the Fourier transform: 

HO = (^«)(0 = 7A3 / < " ' «•:•<• :•</•''• 

For «,»(; 6>(R 3 ) and e > 0, we define the Wigner transform 

s s 

W £ (u,v)(x,£) = (T y ^){u(x+ -y)v(x- -y)), 
W £ (u) = W £ (u,u). 
In the sequel, we denote W £ = W £ (u £ ). 

If (u £ ) is a bounded sequence in Lf oc (M. d ), it turns out that (see [5]), up to 
extracting a subsequence, the sequence (W £ (u £ )) converges weakly to a positive 
Radon measure /1 on the phase space T*R 3 = R 3 . x R| called Wigner measure (or 
semiclassical measure) associated with (u £ ): 

(2.1) G C C °°(M 6 ), Yim(W £ (u £ )^) = J <p{x, f)d/i 

We recall that these measures can be obtained using pseudodifferential opera- 
tors. The Weyl semiclassical operator a (x, eD x ) (or OpY (a)) is the continuous 
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operator from S{R d ) to S'(R d ) associated with the symbol a G S'(T*R d ) by Weyl 
quantization rule 

(2.2) (a w (x,eD x )u)(x) = -^ [ [ a (^±1, eA f(y)e^Mdy. 

We have the following formula: for u,v e S'(R d ) and a E S{R d x R d ), 

(2.3) (W e (u, v),a) s >,s = {v, a w (x, eD x )u) s >,s, 

where the duality brackets (., .) are semi-linear with respect to the second argument. 
This formula is also valid for u, v lying in other spaces as we will see in Section 3. 

2.2. Besov-like norms. In order to get uniform (in e) bounds on the sequence 
(it e ), we shall use the following Besov-like norms, introduced by Agmon and Horman- 
der £Q: for u, f € Lf oc (R 3 ), we denote 




where C{j) denotes the ring {x e R 3 /2 J ' < |ar| < 2 j+1 } for j > and C(-l) is the 
unit ball. 

These norms are adapted to the study of Hclmholtz operators. Indeed, if v is 
the solution to 

—iav + Av +v = f 

where a > 0, then Agmon and Hormander pQ proved that there exists a constant 
C independent of a such that 

\\v\\ B .<C\\f\\ B . 

Perthame and Vega [T2} generalised this result to Hclmholtz equations with general 
indices of refraction. 

We denote for x € K 3 , |x| = v/EjLi x i and ( x ) = ( l + \x\ 2 ) 1/2 - 
For all <5 > |, we have 

(2.4) \\u\\ Lls :=\\{x)- s u\\ L 2<C(5)\\u\\ B ,. 

We end this section by stating two properties of these spaces that will be useful 
for our purpose (the reader can find the proofs in pQ). The first proposition states 
that, in some sense, we can define the trace of a function in B on a linear manifold 
of codimension 1. 

Proposition 2.1. There exists a constant C such that for all f S B, we have 

\ 11/(^1, -)IU 2 (R 2 ) d:c l < CII/lls- 
JR 

The second property gives the stability of the space B by change of variables in 
Fourier space. 
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Proposition 2.2. Let Q%, O2 be two open sets in R 3 , ip : fii — * O2 a, C 2 diffeo- 
morphism, \ G C,i(R 3 ). For all u G B, we denote 

Tu = F~ l (x(u VO)- 

Then 

\\Tu\\ B <C\\x\\clU\\ct\\u\\B- 

2.3. Assumptions. We are now ready to state our assumptions. Our first assump- 
tion, borrowed from concerns the regularizing parameter a £ > 0. 

(HI) a e > e 7 for some 7 > 0. 
This assumption is technical and is used to get a radiation condition at infinity in 
the limit e — ► 0. Next, in order to get uniform bounds on u e , we assume that the 
source terms So and Si belong to the natural Besov space that is needed to actually 
solve the Helmholtz equation 

(H2) HSoIIbJSiIIb <oo. 
It turns out that, in order to compute the limit of the energy source, we shall need 
the stronger assumption 

(H3) (x) N S € L 2 (R 3 ) and (x) N Si € i 2 (M 3 ) for some N > \ + 



3. Bounds on solutions to Helmholtz equations 

In this section, we first establish uniform bounds on the sequence (u 6 ) that will 
imply estimates on the sequence of Wigner transforms (W £ ). It turns out that we 
shall also need to compute the limit of the rescaled solutions Wq and wf defined 
below in order to obtain the energy source in the equation satisfied by the Wigner 
measure fi. 

Before stating our two results, let us define these rescaled solutions. Following 
and 0], we denote 

(3 1) { w °( x ) = e ^ u£ ( ea; )' 

I w\(x) = e^~u £ (qi + ex). 

They respectively satisfy 

f -ia s EWQ + Awq + = S (x) +Si(x - ^-), 
\ -ia e ew\ + Awf + wf = S (x + ^) + Si(x). 

We are ready to state our results on u e , Wq and w\ . 

Proposition 3.1. Assume Sq, S\ G B. Then, the solution u 6 to the Helmholtz 
equation satisfies the following bound 

\\u e \\ B * < C{\\S \\ B + \\SxWb), 
where C is a constant independent ofs. 

Proposition 3.2. Let Wq and w\ be the rescaled solutions defined by Then, 
the sequences (Wq) and (iff) are uniformly bounded in B* and they converge weakly- 
* in B* to the outgoing solutions Wq and w± to the following Helmholtz equations 

f Aw a +w = Sq 
\ Awi + wi = Si, 
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i.e. wo and wi are given in Fourier space by 



Remark. The Helmholtz equation Aw + w — S does not uniquely specify the 
solution w. An extra condition is necessary, for instance the Sommerfeld radiation 
condition. When the refraction index is constant equal to 1, this condition writes 



Such a solution is called an outgoing solution. 

Alternatively, still assuming that the refraction index is constant, the outgoing 
solution to the Helmholtz equation may be defined as the weak limit w of the 
sequence (w s ) such that 



We point out that the two points of views are equivalent in the case of a constant 
index of refraction (which is not true for a general index of refraction). 

We prove the two propositions in the following two sections. As we will see in the 
proofs, our main difficulties are linked to the rays that are emitted by the source at 
towards the point q\ (and conversely). Hopefully, the interaction between those 
rays is "destructive" and not constructive. 

3.1. Proof of Proposition I3TT1 In the sequel, C will denote any constant inde- 
pendent of e. 
The scaling invariance 



First, we note that the bound ||iWoIIb* — C||So||b is established in Agmon-Horman- 
der pQ (see also Perthame-Vega [12] )■ Hence, the proof of Proposition ITTl reduces 
to the proof of the following lemma. 

Lemma 3.3. If a e is the solution to 



(3.2) 




i5w s + Aw s +w s = S(x). 



IKIU- < Kiln-, 

makes it sufficient to prove bounds on Wq. Since Wq is a solution to 

— ia e ewQ + Awq + Wq = Sq(x) + Si(x — -^) 
we may decompose Wq = w^ + a £ , where u>q and a £ satisfy 





-ia £ ea E + Aa E + a e = S x (x - — ) 
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Proof. We want to prove that 

VveB, \(a s ,v)\ <C||Si|| B H B . 
Using Parseval's equality, we write 



(3.3) (a E ,v)= - 

To estimate this integral, we shall distinguish the values of £ close to or far from 
two critical sets: the sphere {|£| 2 = 1} (the set where the denominator in (|3.3|) 
vanishes when e — > 0) and the line {£ collinear to xq} (the set where we cannot 
apply directly the stationary phase theorem to i|3.3|)'l. 

More precisely, we first take a small parameter 5 G]0, 1[, and we distinguish in 
the integral (|3.3|l . the contributions due to the values of £ such that |£ 2 — 1| > S 
or |£ 2 — 1| < 6. Let \ € C^?°(M) be a truncation function such that x(A) = for 
|A| > 1. We denote X s(0 = X^-^r 1 ) ■ We accordingly decompose 

= r + ir. 

First, since the denominator is not singular on the support of \s, we easily bound 
the first part with the L 2 norms 

and using B <— > L 2 , we obtain the desired bound 

(3-4) |/ £ | <C\\Si\\ B \\v\\ B . 

Let us now study the second part II s where the denominator is singular. Up to 
a rotation, we may assume q\ — |3i|ei, where e\ is the first vector of the canonical 
base. We make the polar change of variables 

r sin 6 cos ip 
r sin 9 sin tp . 
r cos 

Remark. In order to make the calculations easier, we write this paper in di- 
mension equal to 3, but the proof would be similar in any dimension d > 3. 

Hence, qi ■ £ = \qi\r sin 9 cost/?, and we get 

C p — i r sin 9 cos ip _ 

II s = / — — ; (Sx$(l-X6))(f;(r,e,<p)y S in6drd0d<p 

J —r z + 1 — iea e v ' 

Now, we distinguish the contributions to the integral d6d<f> linked to the values close 
to, or far from, the critical direction {6 = ^ , <p — 0} (which corresponds to the case 
{£ collinear to qi})- To that purpose, let r\ > be a small parameter and denote 



K = Ur,0, 
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Then if is a compact set. Let k G be such that (1 — xs)k(0, ip) is a localization 
function on if. We write 



ir = 



+ 



-r sin tf cos (p 



—r 2 + 1 — iea L 



-(Siv)(£(r, 0, <p))(l - xs(r))k(0, <p)r 2 sm. 6 drdO dip 



1 — iea.1 



-(Si«)(£(r, 0, <p))(l - x*(r))(l - fc(0, <f))r 2 smOdrd0d<p 



ir = nr + iv E . 

To estimate the contribution III 5 , we apply the stationary phase method. We 
denote a = — \. The phase function is g r {a, ip) = r cos a cos tp so 

dg r 



a 

'9r 

dp 



-r sm a cos ip 



= — r cos a simp 



and the Hessian at the point (0, 0) is 

dV(o,o) = 



at (a, ip) = (0,0), 
at (a,ip) = (0,0), 



-r 
-r 



which is invertible at any point in if. Since If is a compact set, we can apply the 
Morse lemma: there exists a finite covering (J2j)j=i,n ( n G N) of if such that on 
each set Slj, there exists a C°° change of variables (a, ip) >-> (ay, <£j) such that 



2 2 

g r {a,(p) = r-r— ~r—. 



Moreover, we can write (1 — xs)k — Y^j=i Xj where \j G and supp(xj) C fij. 
Then, we make the changes of variables a'j = y/^otj, ip'j — \f\ L Pj- Finally, we 



decompose Xj = XjXj- Thus, we obtain, for the contribution IIP , the formula 



(3.5) 
where 



IIP = £ 



3 i^(-r+af+ V ; 2 ) 

— r + 1 + iea E 



T/5i(r, a'^tp'^Tfvir, a'^ip'^drda^d^. 



T-Sx := r((xjSi) o^(r,a(aj,^),^(ai,<^i))), 



3> 



J" 



1 + Z£ " £ (x 2 w) o £(r, a(aj , ^ ) , p(oy , <^-)) 
— r z + 1 — «ea e J 



d£ 


d(a, <p) 




d(r, a, (/?) 


d(aj,(pj) 


) 



As a first step, using Proposition ^. 21 we directly get T^- Si G £? with 

WTlSiWs < C\\Si\\ B . 
As a second step, we study Tfv. Since for r close to 1, 



-r + 1 + iea £ 



1 — zea F 



< 1, 



we recover, from Proposition ^. 21 



Tfv £ B 



and 



|I?«||fl < C|M|; 
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Now, we apply Parseval's equality with respect to the r variable in the formula (|3.5() 



IIP = Y : T/ 1 5 1 (.-^,.,.)T2«drda'V 7 

^ J -r + 1 + tea £ 3 e ' 3 3 3 

^^^'l^oje-^-^^Sit. - H .))(P - t, ol^'i) 



xJ 7 r ^p(Tjv)(p, a'j , (pj)dtdpda'jdip'j . 

where l{ p >o} denotes the characteristic function of the set {p > 0}. 
Hence, we obtain 



f^/ \\T r ^ p {TjS X (. 



\ur\ > J / .„,r/,s' l (.-i^i J . ) .))(p))|U.dp 



x / ||^(T»(p)|| i2 ^ 



|J// e | < E(/ll^i(p- M)|| L2 dp)( /"||T/^)|| ia dp) 

|i// e | < ^(/ll T /^(/°)IU 2 ^)(/ll T >WIU 2 ^) 



3=1 

|m e | < CHSiHflllwIlfl, 
which is the desired estimate. 

We are left with the part IV s , which corresponds to the directions £ that are 
not collincar to qi . We denote K' the support of (1 — x£){\ ~ k) which is a compact 
set. In K' , we can choose as new independent variables 

»?i = -?i-£, m = |£| 2 - i. 

More precisely, since 

d(77i, 77 2 ) _ ( -qi 



d£ V 2 ^ 

is of maximal rank 2, there exists a finite covering (Clj)j=i, m (m G N) of K 1 such 
that in Clj, we can make the change of variables £ i— > 77. As before, we denote 
Xj = XjXj some localization functions on fl'j such that (1 — Xa)(l — fc) = 2j=i Xj- 
Thus, for j = 1, ... , to, 



-|£| 2 + l + iea% 3 J -rji + iea 
If we denote 



di 
dry 
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and if T\ denotes the Fourier transform with respect to the r\\ variable, ParsevaPs 
equality with respect to r\\ gives 



1 + isa. 



■Sivx'M 



x {T- 1 (ifv ) ) [x ! ) e im / £ dtdx i dr ]2 d m 



Summing over j, we obtain 

\IV £ \ <C\\Si\\ B \\v\\ B , 
which ends the proof of the bound 

\(a £ ,v)\ <C||5i||fl|H| B . 



□ 



3.2. Proof of Proposition I3T21 We prove the result for the sequence (wq), the 
convergence of the sequence (w £ ) can be obtained similarly. As we did in the proof 
of Proposition l3.il we write Wq = Wq + a £ . Since Wq is the solution to a Hclmholtz 
equation with constant index of refraction and fixed source, it converges weakly-* 
to the outgoing solution wo to Aw + w = So. Hence, it suffices to show the following 
result. 



Lemma 3.4. If a e S B* is the solution to 

—ia e ea £ + Aa £ + a £ = Si (x — —) 

x e ' 

then a £ — ► in B* . 

Proof. The proof of this result requires two steps (using a density argument): 

(1) for v e B, we have the bound |(a e ,u)| < C||5i||s||i;||b 

(2) if Si and v are smooth, then (a £ , v) — > 0. 

The first point is exactly the result in Lemma 13.31 It remains to prove the conver- 
gence in the smooth case (the second point above). 
We write 

, e , r e-^<si(om 

{a ,v) = / ; d£_. 

We are thus left with the study of 

(3.6) rm = i g^m < 

where ip = Sxv belongs to 5(R 3 ). 

As in the proof of Lemma 13.31 we distinguish the contributions of various values 
of £. We shall use exactly the same partition, according to the values of £ close to, 
or far from, the sphere |£| = 1 and collinear or not to qi. We shall use the same 
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notations for the various truncation functions. 

We first separate the contributions of £ such that |£ 2 — 1| < 6 and | ^ 2 — 1 1 > S 

e-*^-ty(0(l-X*(0) 



-\£\ 2 + l-isa e 

= r + ir. 



-\£\ 2 + l-iea £ 



In the support of \$, since the denominator is not singular, we can apply the 
non stationary phase method. 
Since qi ^ 0, we may assume q\ ^ and we have 



iq\ 



1 — iea E 



1 — iea £ 



Hence, we obtain the bound 



Since (x*P) and ^lxi' belongs to S, we have, as s — > 0, 

I s -> 0. 

Let us now study the second term II £ . We use the same changes of variables as 
in Section im It leads to the following formula 



3=1 



e -i^(r- a f-v>f) 



-r + 1 + iea e 



Xj (r, a'j , v'M{r, a', , if' Mr da', dip) , 



where 



Xj {r, a j ,(pj) = Xj ° £{r, «(a ■ ,¥?•), ¥>( a j ,<Pj)) —, — r~- i : \ 

r(— r z + 1 — iea e j 



d(a,ip) 



are still smooth functions that are bounded independently from e. 

Using Parseval's inequality with respect to the variables (a), ip'j) for each integral, 

we obtain the bound 



\ir\<ceJ2 

j=l 



—r + 1 + iea £ 



^Xj.niXj^drdXjd^j 



To obtain the convergence of IP , it remains to study an integral of the following 
type 

e s 'w(r) 
; ar, where w € o. 

|r-l|<<5 — ^ + 1 + «£a e 



This is done in the following lemma. 
Lemma 3.5. Vw € S, V6* £ (0, 1), we have 

-flSlir / \ 

e = w\r) 



r\<5 



-r + 



-dr 



iirw(0) + O e ^o(s- & ). 
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Using this lemma, we readily get the estimate 

(3.7) \nr\ < Ce 1 - 9 ye e (0,1), 

which proves that IIP — > as e — > 0. 

There remains to give the 
Proof of Lemma \3.5l We write 



8 _ii2U r / \ 

e e u>(r) 



-<5 



-r + 2ea £ 



dr 



8 -ii2Ji r / \ 

e s w(r) 
/_ 5 r 2 + (ea E ) 2 

= — tea. 



(r — iea e )dr 



~w(r) 



_ s r 2 + {ea e ) 2 



dr - 



rw(r) 



-8 



r 2 + {ea £ } 



- dr 



I + 11. 



We have 



— e - l ^ a - y w(ea e y) 
s y 2 + l 



dy — > — iirw(0), 



and 



77 



(e-^ r w{r) - w(0)) -^--L-^dr + / w ( ) 



r 2 + (ea e ) ; 



'-5 r2 + (ftte) 2 

The last term vanishes because the integrand is odd. Moreover, using the smooth- 
ness of w, we easily obtain that for all 9 £ (0, 1), 



Thus, 



\e-^ r w(r) -w(0)\ < Cg(£j 
"w(r) — w(0)) 



r 2 + (ea e ) 5 



■rdr 



< 



C 



l dr 



and the result is proved. □ 



We are left with the study of IV £ . We use the same change of variables as in 
Section 13.11 



3=1 

m 

= E 

3=1 



-rji + isa £ 



d£ 
dr/ 



ie E 

i=i 



-a, 



,x(W>Xi)(£(r?)) 



d£ 

d.77 



dry. 



The integral obviously converges with respect to all the variables except 771 . It re- 
mains to prove the convergence with respect to the 771 variable, i.e. the convergence 
of 



where 



dm 

-rji + iea £ 



d£ 
dr) 
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is smooth and compactly supported with respect to rj. It is a consequence of the 
fact that the distribution (x + iO)^ 1 is well-defined on K by 

= v.p.(-) — iirS(x). 



x + iO x' 

We conclude that IV s — ► and (a 6 , v) — > as e — * 0. □ 



4. Transport equation and radiation condition on fj, 

In this section, we state and prove our results on the Wigner measure associated 
with (u £ ). Since we established the uniform bounds on (u £ ) and the convergence 
of (u>q), (wf), these results now essentially follows from the results proved in 
We first prove bounds on the sequence of Wigner transforms (W £ ) that allow us to 
define a Wigner measure fi associated to (u E ). Then, we get the transport equation 
satisfied by fi together with the radiation condition at infinity, which uniquely 
determines \i. 

4.1. Results. 

Theorem 4.1. Let So, S\ £ B and A > 0. The sequence (W E ) is bounded in 
the Banach space and up to extracting a subsequence, it converges weak-* to a 
positive and locally bounded measure fi such that 

(4.1) suP^/ / »(x,Odx<%<C(\\So\\B + II^Hb) 2 . 

R>0 n J\x\<RJ^eR 3 

The Banach space is defined as the dual space of the set X\ of functions <p(x,£) 
such that tp(x,y) := J-£—, y (<p(x, £)) satisfies 

(4.2) f sup(l + \x\ + \y\) 1+X \<p(x,y)\dy < oc. 

JR d xGR d 

Theorem 4.2. Assume (HI), (H2), (H3). Then the Wigner measure \x asssociated 
with (u e ) satisfies the following transport equation 

(4.3) £ • V*M = j^(s(x)\So(0\ 2 + S(x - qi )\S\(0\ 2 ) <5(|£| 2 - 1) := Q(x). 

Moreover, fi is the outgoing solution to the equation \4-3\j in the following sense: 
for all test function R £ C£°(R 6 ), if we denote g(x, £) = J Q R(x — £i, £)dt, then 

(4.4) / R(x,£)dti(x,0=- I Q(x,09(x,0dxd£. 
Remark. Here the support of the test function R contains 0, contrary to 

4.2. Proof of Theorem 14. 11 This theorem, that is proved in [2], is a consequence 
of the uniform estimate on the sequence (u E ) in the space B* obtained in Proposition 
13.11 We observe that for any A > 0, 



(4.5) 



||(a : )-i-V( a: )|| ia <C||u e || fl .<C||/|| B) 
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hence, for any function tp satisfying (|4.2|l . we have 
\{W £ {u £ )^)\ _ 

/ [ \ uE \( x + jy)\ uE \( x - fjO / , e a+o/ e \i+oi 1/ u ^ 

< / = f — (x + -y) 2 ^ (x — -y) 2 ^ \tp\(x,y)dxdy 

<C||/||| / sup<M + | y |) 1+( V(x,y)|dy. 

So (W £ (u £ )) is bounded in X^, A > 0. We deduce that, up to extracting a subse- 
quence, (W £ (u £ )) converges weak-* to a nonnegative measure /i satisfying 

(4.6) |<M,<£>| < / sup(|*| + M) 1+ V^)Mz/- 

We refer for instance to Lions, Paul |Hj for the proof of the nonnegativity of /i. 
The bound 1|4.1JI is obtained using the following family of functions 

rf&'V) = -^j2 e ~ MVfl ^({x) < R) 
and letting fj, — > 0, R — > oo. □ 



4.3. Proof of the transport eauation l4.31 This section is devoted to the proof of 
the transport equation satisfied by [i. We first write the transport equation satisfied 
by W e in a dual form. Then, we study the convergence of the source term (the 
convergence of the other terms is obvious). Finally, choosing an appropriate test 
function in the limiting process, we get the radiation condition at infinity satisfied 
by fi. Proving first a localization property, we improve the radiation condition 
proved in 

4.3.1. Transport equation satisfied by W £ . W 6 satisfies the following equation 

is 

(4.7) a e W £ + £ ■ V x W e = -Im W £ (f £ , u £ ) := Q £ . 
This equation can be obtained writing first the equation satisfied by 

v e (x,y) = u £ (x+^y)W(x- |y). 

From the equality 

V y • V x v £ = | [Au £ {x + e -yW(x - £ -y) - AUF(x - £ -y)u £ {x + £ -y)} , 
we deduce 

1 S E 

a £ v £ + iS7 y ■ V X V £ + —[n 2 {x+-y)-n 2 {x~ -y)]v £ = cr E (x,y), 

where 

a E (x, y) := |- [S £ (x + £ -y)W(x - £ -y) - W{x - |y)u £ (x + |y)] . 

After a Fourier transform, we obtain the equation (|4.7|l . 

Then we write the dual form of this equation. Let tp £ S(M. 6 ), we have 

(4.8) a E (W £ , tP) - (W £ ^ ■ V x tP) = (Q £ , tf,). 
By the definition of the Wigner measure fj,, we get 

a E (W £ ,^)^0 and (W £ , £ • V^) -> ( M , £ ■ V*# 
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Hence we are left with the study of the source term (Q e ,ip). 

4.3.2. Convergence of the source term. In order to compute the limit of the source 
term in 1)4. 7JI . we develop 

(Q e ,^) = |Zm((VP(S W),V) + {W s (Sf,U s ),^) 
Thus, the result is contained in the following proposition. 

Proposition 4.3. The sequences (eW e (S e ,u e )) and (eW £ (Sf , u £ )) are bounded 
S'(R 6 ) and for all i/j g S(R 6 ), we have 

(4.9) ]ime(W%S e ,u s ),Tp) s >,s = [ §*(O^(0lKlUK> 

e^O (2tt) j J r3 

(4.10) ]ims{W s (Si,u s ),i;) s ,,s = f ^(O^tO^O^, 
where Wq and W\ are defined in Proposition 1 3. 'A 

Using Proposition 14.31 we readily get 
i 

— \«b , r/ 2 ^ 
1 



lim(Q^) = 57^33^ (y ^(0^0(0^(0,0^+ / Gi(€)Si (0^(9i. O^e 

1^(01^(^-1)^(0,0^ 

|5T(0I 2 ^ 2 - 1)^(91,0^], 



(4tt) 2 



which is the result in Theorem l4.2l □ 
Let us now prove Proposition 14.31 

Proof of Proposition ^. 3\ The two terms to study being of the same type, we only 
consider the first one in our proof. Let tp e S(T*M. d ) and ip(x,y) = ^^(^(x, 0), 
then we have 

e{W e {Sl,u E ),$)s',s = e j Sq(x + ^y)W(x - ^y)<p(x,y)dxdy 

y 

S (x)w%(x + y)tp(e(x + -),y)dxdy. 



Hence, using that ip € S(M. 2d ), we get 

(x + y)0 (x) N (y)» 



s(W%Se,u°U) S ',s\ < C I (x) N \S (x)\ lw ° (x \ y J l \\ + N f xk dxdy 



< C\\(x) N S \\ L qv%\\ B . I sup ^tf^ dy 



x) N (y) h 

for any k > and (3 > 1/2, upon using the Cauchy-Schwarz inequality in x. 
Then, we distinguish the cases \x\ < \y\ and \x\ > \y\ : the term stemming from 
the first case gives a contribution which is bounded by C J ^k-p arL d the second 
contribution is bounded by C J -nxjz- So, upon choosing k large enough, we obtain 
that 



\e(W s (S £ ,u £ ),yj) sr . s \ < C\\(x) N S \\ L 4™ e \\ 



B 
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Now, in order to compute the limit l|4.9|l . we write 

e{W E {S^u e )^) = I Sv{x)W {x + y)L(e(x+y^,y} - <p(0,y)\dxdy 



w £ (x)S (x - y)tp(0, y)dxdy 

= Ie + He- 

Reasonning as above, we readily get that lim £ ^o le = 0. For the second term, we 
have 

II e = Jl4(x)(S *( P {0,.))(x)dx 1 

hence, since Wq converges weakly-* in B* , it suffices to prove that So*<p(0, .) belongs 
to B. We denote 4> = <f(0, ■)■ We have, for ^ > 1/2, 



< C\\S *<f>\\ L 2 =CJ (x) \S o *(/>(x)\ 2 dx 

< CU\\ L i [ (xf\S \ 2 *Mx)dx 



where we used the Cauchy-Schwarz inequality. Hence, we get 

\\So*<j>\\B<C\\(x) N S \\ L z f sup f+f^ dy, 

J R 3 xeR 3 {x) N {y) k 

for any k. As before, this integral converges. Thus, we have established that 
So * 4>(0, .) belongs to B, which implies that 



II e -> / S {x)w^(x + y)ip(0,y)dxdy. 

□ 

4.4. Proof of the radiation condition (|4.4D . It remains to prove that (x satisfies 
the weak radiation condition (|4.4() . 

4.4.1. Support of (J,. In order to prove the radiation condition without restriction 
on the test function R (as assumed in we first prove a localization property on 
the Wigner measure /i. This property is well-known when u £ satisfies a Hclmholtz 
equation without source term. It is still valid here thanks to the scaling of S £ . 

Proposition 4.4. Under the hypotheses (HI), (H2), (H3), the Wigner measure \x 
satisfies 

su PP ( M )c {(*,£) eR 6 / |C| 2 = 1}. 

Proof. Let </> <E C C °°(]R 6 ) and <p £ = <p w (x,eD x ). Let us denote H £ = -e 2 A - 1. 
Since u £ satisfies the Helmholtz equation (|1.1|) . we have 

(4.11) ia e eu £ + H £ u £ = e 2 S e . 

Moreover, H £ is a pseudodifferential operator with symbol |£| 2 — 1. By pseudodif- 
ferential calculus, <ffH £ = OpY (</>{%, £)(|£| 2 — 1)) + 0(e) so, using the definition of 
the measure fj,, we get that 

limU £ H £ u £ ,u £ ) = lim(OpY(ct>(x,t)(\t;\ 2 -l))u s ,u s ) 
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Using the equation (|4. lip , we write 

(<t) £ H £ u £ , u £ ) = e 2 {(j) £ S £ , u £ ) - ia e e(<j> £ u £ , u £ ) = e 2 (W £ (S £ , u £ ),4>) - ia e e(c/) £ u £ , u £ ). 

On the first hand, Proposition 14.31 gives that lim e ^o £ 2 (W £ (S £ , u £ ), </>) = 0. On the 
other hand, (4> £ u £ ,u £ ) is bounded so lim e ^o a e e((f> £ u £ , u £ ) = 0. Therefore, for any 
4> e C C °°(R 6 ), we have J <H|£| 2 - n 2 {x))dn = 0, so supp(^) C {|£| 2 = 1}. □ 

4.4.2. Proof of the condition \4-$ - Using the previous localization property, in 
order to prove the radiation condition Q4.4|) . one may only use test functions R g 
C~(R 6 ) such that supp(i?) C M 6 \U = 0}. 

Let R be such a test function. We associate with R the solution g £ to 

-a e g e + £- V x g e = R(x,£). 

By duality, we have 

(Q £ ,g £ ) = (W £ ,R), 
so that it sufffices to establish the following two convergences: 

(4.12) hm(Q g ,<f) = (Q,g), 

e— >0 

(4.13) lhn{W e ,R) = (f,R), 

e—fO 

where Q et g are defined in Thcorcm l4.2l 

As before, since R € X\ for any A > 0, the limit (|4.13(l follows from the weak-* 
convergence of W £ in 
On the other hand, 

(4.14<jQ e ,<f) = Xm f S^(x)w £ {x + y)g £ {e[x + h,y)dxdy 

+Xm / S 1 {x)w £ 1 (x + y)g £ {q 1 + e[x + ^-},y)dxdy, 
Jr 6 z 

so (Q £ , g £ ) is the sum of two terms of the same type. Such a term has been studied 
in where the following result is proved. 

Proposition 4.5. Assume (w £ ) is bounded in B* and that (w £ ) converges weakly-* 
in B* to w Q . Assume S satisfy (H3). Let R g C^°(M 6 ) be such that supp(R) C 
M 6 \{£ = 0}. Let g £ be the solution to 

-a £ g £ + e-V,.g £ = R(x,t) 

and g(x, £) = J °° R(x + t£, £)dt. Then, we have 

. . ,.. „,.. ,„ _.. , h,y)dxdy = — 

Jr 6 * (27T) J 

Proof. The proof of this result is written in the appendix. □ 
Using the proposition above together with Proposition 13. 21 we get that 
\im(Q £ ,g £ ) 

= i m (-l- f ^(^(^(0,0^ + 7^3/ sliO^iOg^O^ 

V ( 27r ) Jr 3 ( 2 *T Jr 3 



lim / S {x)w £ (x + y)g £ (e[x+^},y)dxdy= j—^ I S (£)^(£)<?(0, £)d£. 



1 



(4tt) 2 



\s m 2 s(e-i)g(o,odc+ \si(o\'s(t- 1)5(91, m 
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Thus, the radiation condition (|4.4(l is proved. 

Appendix A. Proof of Proposition 14.51 
In the sequel, we denote 

G £ = lim / S (x)wq(x + y)g £ (s[x + -],y)dxdy. 

A.l. Bounds on G £ . In order to study G £ , we need a preliminary result on the 
test function g £ . 

Lemma A.l. Let R £ C^°(R e \ {£ = 0}). We denote g £ the solution to 

(A.l) -a e gF+{.V xS f = R{x,t). 

It is given by the explicit formula 

f°° 1 £ 

(A.2) g £ (x, 0=-J exp(-a e |C|~ 1 S )^i?(x - ±s, £)ds. 

Then we have the estimate 

,(x) M Aa- M 



(A.3) VM>0, \gs(x,y)\<C- 



(y) 



M 



where . A . denotes the infimum of two numbers, and C is a constant depending on 
M and R. 

Proof. Let a be a multiindex such that \a\ < M. We denote u> = -Ar. We write 
y a g £ (x,y) = ^(J (^ C ) Q e-^^'^-Rix - us,® da) 

r r+oc 

= dee-*'*/ ds V C {aMeJ) e-^~ s ( l 9 c ) b (-a £ |er 1 s) 

jR3 Js=0 b,c,d,e,f 

(a.4) x (-a^i-hy^n-s^iid^iid^r (J*r) (x - ^ o. 

Using that 

• Re C °°(]R 6 \{|£| =0}) so 

• there exists r , A, B > such that supp(R) C {\x\ < r ] x {A < |£| < B}, 

• R and its derivatives belong to L 1 (M 6 ). 

. i^) b (-a £ icr l S )i<c s , v iei > a. 

• \(id ( ) c (-soj)\ < Cs, V|£|>A 

• in the integral above, s £ [\x\ — r , |x| + r ], so for |x| large enough, we can use 
the equivalence s ~ \x\ ~ (x) (where we dente for a, b > 0, a ~ b if 3 ci, C2 > 
0/ C\a < b < C2tt), we get in (|A.4|I . 



\y a 9 £ (x iy )\ <C(x) M exp(-^(x)). 
The desired estimate follows. □ 
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Using this lemma, we estimate 



\G £ \ < 



/ Wo{x + y)S (x)g e (s(x+ -,y)dxdy 

JR6 1 



<- of Mfe±^ (l + rti^iftwi w- + % A ° r " 

Jr6 (a; + y) 2 + u w 

< cK|| B *||(x)^s (*)IU 2 

/* , \ 1 +0 1 \ n H\x\ + \y\)) M Aa~ M , 
Now, we prove that the integral 



. N (e(\x\ + \y\)) M A a7 M , 
sup (x + y) * +0 (x)- Nl !gj 



/R3 a;GR3 (j/) 

is bounded uniformly with respect to e. 

Let us define the following three subsets in M 6 

a,, A e = {(x,y)GM 6 /M>M}, - {|^l < k 1 " ?/! < 1}, 

( j C £ = {\x\<\y\,\s^ y\>l}, 

where e 1_0 means e 1 ^ 5 with 8 > sufficiently small. 
If (x, y) e A e , then 



(x + y)? +0 (x)- N i (£(|X| + y A a7 - M < C(x)- N ^ +0 ((sx) M A a~ M ) (y)- M . 

Now, we distinguish the relative size of (ex) and e 7 : 

• if (ex) > e~ 7 , we have 

(x)i+°- N ^(ex) M AaJ M ) < (x) l+o-iv l£ - 7 M < £ _ 7 (i +0 -iv 1 +M) i 

• if (ex) < e~ 7 , we have (x) < i(ex), hence we get 

(x)5+°-^((ex) M Aa~ M ) < (x)i +0 -^s^ M < e -7i" e -(7+i)(*+o-"i). 

Now, we choose 1/2 + such that Ni < \ + 0. Then, we get the following bound 
for the contribution of the set A £ to T E 



C£ - 7 M-( 7 +l)(l/2+0-iV 1 ) /" ( y )-M rf2/ _ 



Since TVi > + \ , we can choose M > 3 such that this contribution is uniformly 

bounded with respect to e. 

If {x,y) € B E , then \ey\ < 1 so we obtain 

(x + y)^(x)-^ (£(|X| ± 'f — ^ < C(y)-«+*+°. 

(y) 

Thus, the corresponding contribution to X s is bounded by C J R3 (y)~ M+ ? +0 dy which 

is convergent. 

If (x, y) G C e , then 
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Since \y\ > e 1+0 , if we denote z — ey, we have \z\ > e 5 , for some 5 > 0, which 
implies that 

Thus, we have 

1 E < Ce- 3 [ {-)- M+ ? + %{z) M Ae-< M )dz 



. : ( Z ) M J{z)>el 



. x i +0 (z) M Ae^ M , f , . 1+0j 



where we used the hypothesis (HI) a e > e 7 . 

We distinguish two cases, according to the relative size of (z) and e -7 . We have 
, )i+0 {z) M Ae^ M dz ^ s _ My r (z) k+o-M dz 



and 



(z)<ei V*/ J (z)<e~< 

Hence, we get 

J £ < C i £ .M(l-«)-4-7(|+0) < 

To conclude, we choose M > , which gives that X e is uniformly bounded with 
respect to e. 

A. 2. Convergence of G £ . Wc decompose G £ in the following way 

G e = f w s o (x + y)S (x)(gz(e(x+l),y)-fi(0,y))dxdy 

./B6 i 

Wo(x + y)S'o(a:)((r e (0, y) - g(0, y))dxdy 

w o( x + y)So{x)g(0, y)dxdy 
= I e + II e + IJI e 

Using the same method as in Section TA. II we prove that I £ , II £ — ► 0. Then, we 
may write 

III e = [ w £ o (x)(S o *g(0,.))(x)dx. 

JlSL e 

Moreover, we established in the proof of Proposition ^. 3l that if 4> is rapidly decreas- 
ing at infinity, then So * (f> G B. Hence, since (wq) converges weakly-* in B* to wo, 
we get 

lim IIL = [ w (x)(S * 5(0, .))(x)dx, 

which ends the proof. 
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